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Abstract. A wide class of non-autonomous nonlinear parabolic partial differential 
^ ! equations with delay is studied. We allow in our investigations different types of delays 
^ I such as constant, time-dependent, state-dependent (both discrete and distributed) to be 
presented simultaneously. The main difficulties arise due to the presence of discrete state- 
dependent delays since the nonlinear delay term is not Lipschitz on the space of continuous 
functions. We find conditions for the local existence, uniqueness and study the invariance 
. principle. 

< 

^ ; 1 Introduction 

We consider non-autonomous parabolic partial differential equations (PDEs) with delay. 
Studying of this type of equations is based on 1) the well-developed theories of the delayed 
ordinary differential equations (ODEs) [TOl[IIlE] and 2) PDEs without delays [H El [H US] . 
Under certain assumptions both types of equations describe a kind of dynamical systems 
CN| ■ that are infinite-dimensional, see [U [2S1 E] and references therein; see also [301 El IISI El II] 

and the monograph [M] that are close to our work. 
^ . In evolution systems arising in applications the presented delays are frequently state- 

! dependent (SDDs). The theory of such equations, especially the ODEs, is rapidly devel- 
oping and many deep results have been obtained up to now (see e.g. [20l |3ll |32l [17] and 
also the survey paper [12] for details and references). The PDEs with state-dependent 
delays were first studied in [2T1 [T3l 122] . An alternative approach to the PDEs with discrete 
^ . SDDs is proposed in [23]. Approaches to equations with discrete and distributed SDDs 



, are different. Even in the case of ODEs, the discrete SDD essentially complicates the 
study since, in general, the corresponding nonlinearity is not locally Lipschitz continuous 
on open subsets of the space of continuous functions, and familiar results on existence, 
uniqueness, and dependence of solutions on initial data and parameters from, say [HI E] 
fail (see [33] for an example of the non-uniqueness and [l2j for more details). It is impor- 
tant to mention that due to the discrete SDDs such equations are inherently nonlinear. 
In this work, in contrast to previous investigations, we consider a model where two dif- 
ferent types of SDDs (discrete and distributed) are presented simultaneously (by Stieltjes 
integral). Moreover, all the assumptions on the delay (see (Al)-(A5) below) allow the 
dynamics when along a solution the number and values of discrete SDDs may change, the 
whole discrete and/or distributed delays may vanish, disappear and appear again. This 
property makes it possible to study "flexible" models where some subsets of the phase 
space are described by equations with purely discrete SDDs, and others by equations with 
purely distributed SDDs, and there are subsets which need the general (combined) type 
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of the delay. A solution could be in different subsets at different time moments. This 
property particularly means that not only the values of the delays are state-dependent, 
but the type of the delay is state- dependent as well. 

The first goal of the present paper is to study the basic properties of solutions - the 
existence and uniqueness as well as to extend the fundamental invariance principle to the 
case of PDEs with discrete SDDs. The second goal is to attract attention of researchers 
from such fields as, for example, mathematical biology and physics to this wide class of 
delay equations and emphasize that the crucial assumption on the delay (see (A5) below) 
is an "inner property" of the delay which could be successfully used in a wide range of 
other delay systems. We hope that our results provide a basement for further study of 
qualitative (asymptotic) properties of solutions. 

The existence and uniqueness results for a particular case of autonomous systems 
were announced in |25] . For a survey of the existing literature on the invariance principle 
see [26]. In the present paper the emphasis is on the delayed term, not on the partial 
differential operator. To the best of our knowledge the invariance principle for PDEs with 
SDDs has not been studied before. 



2 Formulation of the model and examples 

Let X be a Banach space with the norm || ■ ||, let r > be a constant. Denote by 
C = C([— r, 0]; X) the space of continuous functions ip : [— r, 0] — )■ X with the supremum 
norm 1 1 ■ 1 1^;. As usually for delay equations [TOl [11], for any real a < 6, t G [a, h] and any 
continuous function u : [a — r,h] — )■ X, we denote by Ut the element of C defined by the 
formula Ut = Ut{9) = u{t + 9) for 9 G [— r, 0]. Consider an infinitesimal generator A of 
a (compact) Cq semigroup {e~^*}t>o = {T{t)}t>Q on X satisfying ||T(t)|| < e"* for all 
t > 0, where a; G M is a fixed constant. 

We are interested in the following non-autonomous parabolic partial differential equa- 
tion with state-dependent delays (SDD) 

^^ + Au{t) = B{t,Ut), t>a (1) 

with the initial condition 

Ua = U\[a-r,a] =^eC = C([-r,0];X). (2) 

The delay term i? : M x C — )■ X has the form 

=G(t,7/;(0),F(t,V^)), (3) 

where G:MxXxX— J-Xisa continuous mapping and the delay functional F : 
M X C — !■ X is presented by a Stieltjes integral (simultaneously includes discrete and 
distributed SDDs) 

F{t, f p {t, ■ dg{9, t,^), p:RxX^X. (4) 



Assumptions on g are formulated below (see (Al)-(A5)). 
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The class of equations described by (JI]),©,© is very wide and includes many equa- 
tions which were intensively studied during past decades. Below we mention just two 
examples and refer the reader to [SU [21] for more references and discussion. 

Example 1. Let X = where Q C is a smooth bounded domain. Operator A 

is a densely-defined self-adjoint positive linear operator with domain D{A) C L^(f2) and 
compact resolvent, which means that A : D{A) — )■ generates an analytic semigroup. 

If we choose G(t,u,v) = v — du {d > is a. constant), then B(t,Ut^ = F{t,Ut) — du and 
equation ([1]) reads as 

d 

—u{t, x) + Au{t, x) + du{t, x) = {F{ut)) (x), (5) 

with, for example, 

{F{^)){x)^ j' yj{^{e,y))fix-y)dyYdgie,^), x e Q, (6) 

where f : Q — Q ^ is a bounded measurable function, p : M — )■ M. This non-local 
autonomous equation is studied in [25j. It is clear that the integral delay term given by 
(E]) includes the cases: 

a) purely discrete SDDs: {F{^p)){x) = Y.k InP ^'^^~'^k{i'),y)) f{,x - y) dy; 

b) purely distributed SDD: {F{i))){x) = jl, {J^p{ip{0,y)) f{x - y)dy}-^{e,ip) dO. These 
cases have been studied in [211 [221 [23] . 

Similarly, one may consider /oca/ delay terms (discrete and/or distributed SDD) 

{F{ij)){x)^ [\{ij{9,x))-dg{e,i:), x e Q. (7) 

The above type of equations includes the diffusive Nicholson's blowflies equation (see e.g. 
[28] ) with state-dependent delays, i.e. equation ([5]) where —A is the Laplace operator 
with Dirichlet or Neumann boundary conditions, Q C is a bounded domain with a 
smooth boundary, the nonlinear (birth) function p is given by p{w) = pi ■ we~^ , pi G M. 

Example 2 (reaction-diffusion system with delay). Suppose Q C is a 
bounded region with a smooth boundary dfl, (9„ is the outward normal derivative on 
dQ, A is the Laplacian operator on Q. Consider the system 

' |f(x,t) = diAu'{x,t) + G'{t,F%t,ul{x, ■),... ,uT{x,-))), t > a, x e n, 
< a'{x)u'{x, i) + dnu\x, t) = 0, t > a, x G dVt, (8) 

^ M^(x, a + e) = ip\x, e), e e [-r, O], x G fi, 

where i = 1, . . . ,m. In di > and = we agree that no boundary condition 
applies to u\ a* G C^~^°' (d^l) , a G (0,1). Functions : — )■ R are locally Lipschitz 
and the delay functionals : R x C{[—r, 0]; R™) — > M are presented by Stieltjes integral 
(simultaneously includes discrete and distributed SDDs) similar to The system 
could be presented in the form ([I])-([2]) as follows. We set X = (7(0; R™) and (see e.g. 
[T9| p. 5], |18] and references therein) let be the operator defined by A^yi = diAyi 
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(or A^^Vi = if = 0) on the domain D{A^) = {yi G C\Q) n C^iQ); aY + dnV' = 
on dn} (or L)(A°) = C(0) if di = 0)). The operator Ai is the closure of on C(0) 
and y4 = (Aj)™^. We denote {Tj(t)}t>o the Co-semigroup on (7(0) generated by Ai and 
T = (Ti)™i. It is well-known [IH] that T is a Cg-semigroup on X = ^(0; R*^) that is 
analytic (and compact if all di > 0) and A is its generator. The system ([H]) was studied 
(without state-dependent delays), for example, in [,18j. 

As an application one can consider the n-species Lotka-Volterra model of competition 
with diffusion and delays given by 



diAu'{x, t) + hiU'it, x) I - YJj=l (^ij I-r "^(^' ^ + ^)) (^9iji(^, Ut) 



X E fl, 



X E dVl, 



u'{x,a + e) = (p'{x,e), 0G[-r,O], 

(9) 

where 6,, Cij are positive constants and gij are nondecreasing with respect to the first 
coordinate and gij{0, ■) — gij{—r, ■) = I. Many interesting properties of this system (au- 
tonomous and without state-dependent delays) were discussed in [19] (see also references 
therein) . 

The approach developed in the present article is applicable to more general classes of 
equations of the form ([T]) with the nonlinearity B, for example, as follows (c.f. (I3l)) 

B{t,ut) =G{t,F\t,Ut),...,F\t,ut)), 

with be as in (jl]). We formulate our results for B given by (13]) for the simplicity of 
presentation and motivated by (Q. 



3 Local existence and uniqueness 

The following assumptions on the time- and siate-dependent delay are generalizations to 
the non- autonomous case of the ones proposed in [25] . 

(AI) For any {t,(f) G M x C, the function [— r, 0] 3 g{-,t,(f) W is of bounded 
variation on [— r, 0]. The variation V^j.g of g is uniformly bounded i.e. 
3Mvg > : V(t, ip)eRxC V^,g{-, t, if) < Mvg. 

It is well-known that any Lebesgue-Stieltjes measure (associated with g) may be split 
into a sum of three measures: discrete, absolutely continuous and singular ones. We will 
denote the corresponding splitting of g as follows 

g{9, t, v?) = gdi9, t, v?) + 9acid, t, (p) + gs{0, t, if) = gd{0, t, v?) + g^O, t, y?), (10) 

where gd{0,t,{p) is a step-function, gac{0,t,ip) is absolutely continuous and gs{0,t,{p) is 
singular continuous as functions of their first coordinates (see [16] for more details) and 
we denote the continuous part by gc = gac + Qs- 
Our next assumptions are 

(A2) For any 9 G [— r, 0], the function g^ is continuous with respect to its second and 
third coordinates i.e. MO G [-r,0], V(t,(^), (r,^^") G M x C : (t",(^") {t,ip) inRx 
C(n^ +oo) g,ie, r, v9") ^ g,ie, t, if). 
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(A3) The step-function gd{6,t,(f) is continuous with respect to {t,(f) in the sense 
that discontinuities of gd{0,t,ip) at points {6k} C [— r, 0] satisfy the property: there are 
continuous functions : M x C — )■ [0, r] and /i/^ : M x C — M such that 6^ = —rjkit^ip) 
and hk(t,ip) is the jump of gd at point 6k = —rik(t,(p) i.e hk(t,ip) = gd{6k + ^,t,ip) — 
gd{6k - 0,t,ip). 

Taking into account that gd may, in general, have infinite (countable) number of points 
of discontinuity {6k}, we assume that the series Ylk^kit^f) converges absolutely and 
uniformly on any bounded subsets of M x C. 

Following notations of ffTOj) . we conclude that (A3) means that for any (t, x) G M x C 
one has $d(t, x) = x(^) dgd{6, t, ip) = xi^k) ■ hk{t, (p) = x{-Vk{t, ^)) ■ hk{t, ip). 
Here all rjk and hk are continuous functions. 

The first result is (c.f. |25l lemma 1]) 

Theorem 1. Assume G : M.x X x X X is a continuous mapping andp : RxX — A 
(see is Lipschitz (| |p(t, u) — p{s, v)\\ < Lp{\s — t\ + — i;||), satisfying \ \p{s, m)| | < 
Ci||m|| + C2,V(s,u) G M X a with Ci > 0. Under assumptions (Al)- (A3), the nonlinear 
mapping i? : M x C — t- A, defined by is continuous. 

Remark. It is important that nonlinear map B is not Lipschitz in the presence of 
discrete SDDs. The last means that discrete delays may be present, but be constant or 
time- dependent only (i.e. gd{6,t,!f) =^d{6,t)). 

Proof of theorem 1. Since a composition of continuous mappings is continuous, it is 
enough (see ([3])) to show the continuity of F defined by (jl]). 

We first split our g in continuous and discontinuous parts gc = gac + gs and gd, 
respectively (see (fTOj) ). This splitting gives the corresponding splitting F = + Fd, 
where F^ corresponds to the continuous part gc = gac + Qs- 

Case 1. Let us first consider the part F^. We write 

Fc{t\<f)-Fc{t\^) = h + l2, (11) 

where we denote 

h = h{if,^)^ j' [p{t\m)-p{t\m)\ dgc{6,t\^), (12) 

h = hi^, ^) = r p{t\ m) d [gc{6, t\^)- gc{6, t\ ij)]. (13) 

J — r 

Using the Lipschitz property of p and (Al), one can check that 

\\h\\<L,{\t' -e\ + \\^-i,\\c)-Mvg. (14) 

This shows that | | — )■ when — + | |y9 — '?/'| |c — > 0. To show that | I/2I | ^ (when 
t^ — )• t^ and —)■'?/; in C) we use assumptions (Al) and (A2) to apply the first Helly's 
theorem [161 page 359]. 

Case 2. Now we prove the continuity of Fd [discrete delays). Let us fix any G C, t G M 
and consider any sequences C C and {t^} C M such that 1 — |c — > and t^ ^ t 
when n — )■ cxD. Our goal is to prove that ||Fd(t"v^'^) — Fdit, ip)\ \ — )■ 0. 
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Following the notations of (A3) we write 



k 

and remind that it could be a series or a finite sum. We split as follows 

F,{r, - F,{t, (^) = Kl' + + G X, (15) 



where 



^ J2 {pit"", ^'\-Vk{t\ y^"))) - Pit, V{-Vk{t\ y.")))} ■ ^^), 

k 

^ J2pit,vi-Vk{t\v"))) ■ [hk{e,vn - h{t,^)] , 

k 

^ ^"))) - V{-Vk{t, ^)))} ■ h,{t, ^). 

k 

Using the Lipschitz property of p one may check that 

lli^ril < L,{r -t\ + \\^^- ^\\c) ■ \hk{t'',^'')\. (16) 

k 

Now we discuss The growth condition of p implies \\p[t,ip{—rik{t^ < 

(Cillv^llc + C2)- Hence 

\\K^\\ < {CilMc + C2) -YMf^'Vl - hk{t,^)\. (17) 

k 

In a similar way we obtain 

\\K^\\ < LpJ2\hkit,v)\ ■ M~Vk{r,vn)-vhVk{t,v)\\. (18) 

k 

Now we show that H-f^'j'^ll — ?■ as n — t- 00 for j = 1,2,3. The first property H-R'"!! — > 
follows from (A3) and f|T6l) . In f|T7|) . the series converges uniformly with respect to n since 
the condition \ \ip'"- — {p\\c + — t| — )■ implies that {{t, (p), (t", v^")} is a bounded subset 
of M X C. Assumption (A3) guarantees that each \hk{t^, </?") — hk{t, (p) \ is continuous with 
respect to (t", (p^) and tends to zero when n — )■ cxd. Due to the uniform convergence of 
the series in f|T7j) (see (A3)), we arrive at H^^^^ || 0. To show that H-ft'sH — )■ we also 
mention that each \hk{t,ip) \ ■ \ \(p{—rik(t"-,ip"')) — (p{—rik(t,ip)\\ (see (fTSj) ) is continuous with 
respect to (t", (p^) and tends to zero as n — 00 due to (A3) and the strong continuity of 
(p E C. The uniform convergence (w.r.t. (t",(y?")) of the series in (fTSl) follows from the 
estimate \hk{t,ip) \ ■ | |v3(-r/fc(r, </?")) - ip{-r]k{t,ip)\\ < \hk{t,ip)\ ■ 2||v?||c (the right-hand 
side is independent of n!) and the Weierstrass dominant (uniform) convergence theorem. 
We conclude that \\K^\\ — )■ 0. Since all \ \Kj'\\ — )■ as — )■ 00 for j = 1, 2, 3 we proved 
the property \ \Fd{t^, v?") — F^it, {p)\ \ — )■ 0. We shown that both Fc and F^ are continuous. 
The proof of theorem 1 is complete. I 

In our study we use the standard 
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Definition 1. A function u G C{[a — r,T]; X) is called amild solution on [a — r,T] 
of the initial value problem (QP, ^ if it satisfies ^ and 

«(t) = e-^(*-"V(0)+ / e-^(*^^)5(s,M,)rfs, te[a,T]. (19) 

J a 

Theorem 2. Under the assumptions of theorem 1, the initial value problem (QP, ([^ 
possesses a mild solution for any E C . 

The existence of a mild solution is a consequence of the continuity of -B : M x C — )■ X, 
given by theorem 1, which gives us the possibility to use the standard method based on 
the Schauder fixed point theorem (see [3, theorem 3.1, p. 4]). 

Theorem 3. Let all the assumptions of theorem 1 are valid. If additionally 
\\G{t,u,v)\\ < ki{t){\\u\\ + \\v\\) + k2{t) with ki are locally integrable on [a, oo), then 
a mild solution is global i.e. defined for all t > a. 

The statement follows from theorem 2 and [M!, theorem 2.3, p. 49]. 

To get the uniqueness of mild solutions we need the following additional assumptions. 
(A4) The total variation of function Qc = Qac + Os satisfies 

3Lvg^> :yt\t' > a ^V^M-,t\ip) - g,{;t',^)] < LvgMt' -t'\ + \\v-^\\c)- (20) 

(A5) The discrete generating function gd satisfies the following uniform condition: 

• there exists continuous function r]ign{t) > 0, such that all rj^ and "ignore" 
values of {p{9) for 9 G {—riign(t),0] i.e. 

^v^9n{t) >o:yt>a,\/ip\ip^eC:yee [-r, -v^9n{t)], ^\e) = ^\e) =^ 

VA; G N ^ r]k{t,ip^) = r]k{t,ip'^) and hk{t,ip^) = hk{t,ip'^). 

Remark. Assumption (A5) is the natural generalization to the non- autonomous case 
of multiple discrete state- dependent delays of the condition introduced in [23]. In f2R 
the function rjignit) was constant rjignit) = rjign > 0. For more details and examples see 
JM^ and also 

Theorem 4. Assume (Al)- (A5) are valid, p is as in theorem 1, mapping G : MxXx 
X X is continuous and locally Lipschitz with respect to its second and third coordinates 
i.e. for any R> there exists Lg,r > such that for all t > a, < R one has 

\\G{t,u\v') - Git,u',v')\\ < Lg,r {\\u'-u'\\ + \\v'-v'\\) . (21) 

Then initial value problem ^ possesses a unique mild solution on an interval of the 
form [a, b) where a < b < +oo for any ip E C. The solution is continuous with respect to 
initial data i.e. Hy?" — ip\\c ^ implies — Ut\\c — ^ for any t G [a, 6). Here is 
the unique solution of (QP, ^ with initial function yj" instead of (p. 

Proof of theorem 4- For the simplicity, we first consider a particular case when the 
generating function g = g^ = gac + gs i-e. F = Fc does not contain the discrete delays. 
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Let (t^,ip), belong to a bounded subset B C M x C. We use the splitting (ITT]) . 

One can see that (see f|T3|) ) 

\M<Mis- V'>Mt\y^) - 9c{tMj)]. (22) 

Assumption (A4) and f|T^ imply that Fc is locally Lipschitz i.e. for any R > there 
exists Lp^ R > such that for all \ \(p\\c < R, Wi^Wc < R, a < < a + R one has 

\mt\<P) - F,(t2,V^)|| < Lp^^n {\t'-t'\ + \\^-i;\\c). (23) 

Consider a sequence {v'"} C C such that ||<y9" — <y9||c — )■ as n — )■ oo. Denote by 
u(t) = u(t; (p) any mild solution of ([1]), ([2]) and by u"'{t) = u^{t; y?") any mild solution of 
([T]), ([2]) with initial data v?" £ C*- The existence of these solutions is proved in theorem 
2. The Schauder fixed point theorem (see e.g. P?i theorem 2.1, p. 46]), used in the proof 
of theorem 2 implies that one can choose i? > and T G [a, &) to have ||v5"||c ^ -R; 
||M"(t;yp")|| < for all t G [a,T]. 

Using the local Lipschitz property of mapping G (pTj) . (!23l) and the form ([3]), we get 

||S(t,<^) - B{t,i;)\\ < Lg,r{1 + - ^llc + LG,ij||F,(t,(^) - F,(t, V>)||. (24) 

Hence for any t G [a, T] one has (we remind that ||T(t)|| = Ije""^*!! < e"^* and F = Fc) 

\\ut-u'l\\c< e"(^-'^) 1 - (^"1 Ic + Lg,r(1 + L,.,,R)e-(^-'^) ■ f \\us - u^Wc ds. 

J a 

The last estimate (by the Gronwall lemma) implies 

\\ut - <||c < e^C'«(i+^^-«)""'^""'e'^(^-'') -W^- </?"||c. 

That is 



\\ut - <||c < Ct ■ ||</^ - </^1|c, Vt G [a,T], Ct = e'^^^"") exp{LG,R(l + L^,,,i)e'^('^-'^)}. 

(25) 

It proves the uniqueness of mild solutions and the continuity with respect to initial data 
in the case g = g^. 

The second particular case g = gd (the purely discrete delay) and only one point of 
discontinuity has been considered in detail in |23] (the autonomous case). It was proved 
in [23] that (A5) implies the desired result. 

Now we consider the general case (both discrete and distributed delays, including the 
case of multiple discrete SD-delays). 

Using the splitting F = Fd + Fc, we have, by definition of mild solutions, 

u^{t) - u{t) = e-^(*-")(<^"(0) - ^(0)) + Te-^^*-^) {i^d(r,<) - F,(r,n.)} dr 

J a 

+ Te-^^*-^) {i^c(r,<) - Fc{T,Ur)} dr. 

J a 

Using fl23l) . one gets for all t G [a, T] 

\\u^{t)-um < ||<^"(0) -<^(0)||e'^(^-") +e-(^-»)LG,ij f \\Fd{ry^ - F,{r,uMdr 

J a 
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J a 

+ LgA^ + LF^,n)e''^^~^'> f \K - UrWcdr, (26) 

J a 



where 



G"(t) = ||^"(0) - <^(0)||e'^(^-'') + LG,ije-(^-'^) T ||F,(r,<) - F,(r,nOII cir (27) 

J a 

is a nondecreasing (in time) function. Multiply the last estimate by e~*'^'^'^'^^+'^^='^)^'^''^ 
to get 

which, after integration from a to t, shows that {G^{t) is nondecreasing) 

J a 



We have 

ft 



LgA^ + ^F.,ii)e"(^~'^) ^ IK - UrWcdr < G"(t) (^e(*-'^)^«.«(i+^^-«)^"*^""' - l) . 

We substitute the last estimate into to obtain 

IK - Utile < ^"(t) ■ e^^-'^^^GM^+LF^.R)^-^''-''\ (28) 

Let us fix any c> b and denote by rjign = m.m{riign(t) : t G [a, c]}. By Assumption 
(A5), rjign{t) > and continuous, so rjign > 0. Let us denote by a = min{6, a + rjign} > a. 

Now our goal is to show that for any fixed t G [a, cr) one has G^{t) — )■ when n — )■ oo 
(we remind that | ly?" — (/9| |c — ^ 0). 

Let us consider the extension functions 



^(0) se(0,a) ^"""^ ^ 



<^"(s) sG[-r,0]; 
<^"(0) se{0,a) 



An important consequence of (A5) is that Fci(t,Ut) = Fd{t,(pt) for all t G [a, a) and 
any solution u : [a — r,a) ^ X , satisfying Ua = (f- In the same way Fdit, m") = Fd{t, ip^) 
for all t G [a, a). Hence the continuity of Fd implies ||Frf(r, ■^") — Fd(r, ^^)|| — )■ for any 
T G [a, a). 

Remark. We notice that the case we consider now is simpler that the one in the proof 
of theorem 1 (see (173]) ) since we estimate Fd at the same first coordinate (time moment 
T G la, a)). 
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The property ||Fd(r, — Frf(r, ^^)|| — )■ for any r G [0,0") and the uniform bound- 
edness of the term allows us to use the classical Lebesgue-Fatou lemma (see |35l p. 32]) for 
the scalar function \ \Fd{T,'ip^) — Frf(r, ^^)|| to conclude that G"'{t) — )■ when tt, — )■ 00 (for 
any fixed t G [a, a)). Hence fl28|) gives the continuity of the mild solutions with respect 
to initial functions for all t G [a, cr). Particularly, it gives the uniqueness of solutions. 
For bigger time values we use the chain rule (by the uniqueness) for steps less than or 
equal to, say (a — a)/2 (for more details see [23]). Since the composition of continuous 
mappings is continuous, the proof of theorem 4 is complete. I 

Remark. Discussing the proof of theorem 4, we see that in the case g = gc (no 
discrete state- dependent delays) the delay mapping F = is locally Lipschitz continuous 
(see ^EM)) ^ consequence, one has a standard estimate for the difference of two 

solutions in terms of the difference of initial functions \\ip — v^^Hc- In case of the 
presence of discrete SDDs we have estimate [2^1 with G" defined by [2l\ ) and for the 
Lebesgue-Fatou lemma it was enough to have property ||Fd(r, ^") — Fd{T,'i^^)\ \ — )■ which 
does not provide information on the difference \ \F{t, ip"') — F{t, (p)\ \ in terms of\\(f — (f"'\\c 
(it is definitely not a Lipschitz property). A way to get such an information is to use the 
modulus of continuity Uf{6;Y). We remind that cOf{6;Y) = sup{||/(x) — f{y)\\ ■ x,y G 
Y, \ \x — y\ \ < 6}. For the simplicity of presentation we consider Fa with one discrete SDD. 
We have (see / I73]) with t = t^ = t) 

F,{t, ^") - F,(r, ^) = p(r, ^"(-r7(r, ^"))) ■ /i(r, ^") - p(r, ^(-r/(r, ^))) ■ /i(r, ^). (29) 

The estimates for Kf, i = 1,2,3 (see [T^) - 17^) with k = 1) show that 

||F,(r,yp") - F,(r,<^)|| < LpV^.g^ ■ ||<^ - <^"||c 

+{CM\c + C2) ■ (||</^-¥^"||c;V') (uj^ ;[-r,0]) , 

where we denoted by Y = {(t, (^t), (t, ^9") : t G [a, cti],??, G N}, with a < ai < a . By (Al) 
one has V^^gd < Myg and 

\\Fd{T, (^") - F,(r, <^) 1 1 < Lj,Mvg ■ [| |v? - V?"l |c + [^r, (\W - (^"l |c; ; [-r, 0]) 

+ (Ci||^||c + C2)-a;/.(||<^-<^"||c;:^) • (30) 

Since (f'^^^pin C, we see that Y is compact. Using (A3) and the classical Cantor 
theorem, we know that h and r] are equicontinuous on Y and Uh (^\\(p — (p"'\\c]Y^ — ?■ 

and (^\\(p — v^"! lei as \\ip — (p"'\\c 0. We remind that Fd{t, Ut) = F^it, ^pt), 

Fd{t,u^) = Fd{t,p^) and use \\pt — '^t'Wc ^ W'^ — V^Wc- Finally, one can substitute the 
estimate ^ into [2l\ ) and then the estimate for into l[28\) to get an estimate for the 
difference of two solutions in terms of the difference of initial functions \ \(p — (f^\\c. 

Remark. In the theory of ordinary differential equations with SDDs it is usual to 
restrict the class of initial functions ip to Lipschitz ones f31\ [7 ^ . In this case one restricts 
the set of SDDs (both rj and h) to Lipschitz mappings. In such a situation the previous 
remark evidently provides the Lipschitz property of Fd and hence F . It follows from the 
property f G Cip{Lf; Y) =^ ^fi.^] Y) ^ Lf ■ ^ and the estimates above. 
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4 Invariance 



This section is devoted to an extension of the fundamental invariance principle [18] to the 
case of PDEs with discrete state-dependent delays. 

We also refer the reader to [26] for important generalizations of the invariance principle 
in several directions (without SDDs) and for a survey of the existing literature on the 
subject. In the present paper the emphasis is on the delayed term, not on the partial 
differential operator. 

Following [18], we assume that the next hypotheses are satisfied: 

(HI) D is a closed subset of [a — r, oo) x X and D{t) = {a; G X : (t, x) G D} is nonempty 
for each t > a — r. 

(H2) V is the closed subset of [a, oo) x C defined hj V = {(t,ip) : 
ip{e) G D{t + 9) for all -r <e< 0}. Also, V{t) = G C : (t, ^) G V} for each 
t > a, and we assume that I^(t) is nonempty for each set t > a. 

(H3) For each b > a there are a K{b) > and a continuous nondecreasing function 
rjb : [0, 6 — a) — )• [0, oo) satisfying //^(O) = with the property that if a < ti < ^2 < 
b,Xi G D(ti), and X2 G D{t2), then there is a continuous function w : [ti,t2] — ^ ^ 
such that w{ti) = Xi,w{t2) = X2,w{t) G D{t) for ti <t < t2, and 

\w{t) - w{s)\ < r^,{\t - s\) + k(b)\t - 

12 — ti 

for all s, t G [^1,^2]- 
(H4) B is continuous from D{B) into X where V C D{B) C [a, 00) x C. 



Remark fT8^, page 16]. If D is convex then (H3) is automatically satisfied by defining 



{t2 - t)xi + (t - ti)x2 . ^ . ^ . 

w{tj = — for ti <t < ^2- 



it2 - tl) 

{s—t)xi + (t—s)x2 
(t2-tl) 



< 



We see that \\w{t) — w{s)\\ - 
hence (H3) is satisfied with K{b) = 1 and r]f, = 0. 
We use the notation 



(t2-tl) 



\t — s\ for tl < s < t < t2 and 



d{x;D{t)) = M{\x-y\:yeD{t)} for xeX,t>a. 

The fundamental criterion for the invariance of the set V, called the subtangential condi- 
tion, (see [IHl (2.2)]) is given by 

lim ^d(e-'^^ip(0)+ 6-'^^'+''-''^ Bit, if) ds;D(t + h)] for (t,ip)eV. (31) 
h-,o+ h \ Jt ) 



The following result is an extension of [TH| theorem 2] to the case of general state- 
dependent delay. 
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Theorem 5. Assume (Al)- (A5), (H1)-(H4) and l[JJ\) are valid. Let mapping p he 
as in theorem 1 and G:MxXxX— J-X satisfy the property: for each R > there are 
an Lg^r > and a continuous vr : [0, oo) — )■ [0, oo) such that z//j(0) = and 

\\G{t\u\v') - G{t^u\v')\\ < UR{\t' -t'\) + Lg,r {\\u' - u'W + \\v'- v'\\) (32) 

for all \ \v^\ \ < R and a <t^ ,t^ < a + R. 

Then initial value problem ([IP, ^ has a unique mild solution on an interval of the 
form [a, h) where a < b < +oo for any ip E C. If additionally ip G T>{a), then u{t) G D{t) 
for a <t < b and if b < +oo then \ \ut\\c — J-oo as t ^ b — 0. 

Remark. In contrast to fTR theorem 2], the nonlinear term B in equation (Qp is not 
Lipschitz (with respect to the second coordinate, c.f. fl^ property (2.3), p. 18]) in the 
presence of discrete SDD. So fTR theorem 2] could not be applied to our case. Instead, 
assumption (A5) provides the uniqueness of mild solutions and saves the line of the proof 
presented in 118]. 

The proof of theorem 5 follows closely that of theorem 2]. The last consists of 
eight lemmas and the final part which spends pages 35-43 of the original article. There 
is no need to repeat these lemmas since they are not affected by the lack of the Lipschitz 
property of B and we refer the reader to [18] for all notations and details. Here we only 
remind the main steps of the original proof and give the new part of the proof based on 
assumption (A5). 

First, for fixed a > a,eo > and any e G [0,£o] the e-approximate solution w is 
constructed. It is done by a careful construction (see [IB] for all details) of an increasing 
sequence {tj}^ C [a, a + eo] such that w{ti) G D{ti) and (see [HI (4.6)]) 



<e{t,+^-ti). (33) 



Let be a decreasing sequence such that — )■ as n — i- oo and for each n > 1 

let and {t"}^Q be as constructed above with e = en,ti = t", and w = w". Denote by 
7" : [a, a] — )■ [a, a] the function 7"(t) = t" whenever t G [^",^^+1]- 

Remark A. In addition to consideration in fW] . we assume that a + Eq — a < rjign = 
m.mt^[a,c] vit) /o?^ some fixed c > a. Since we prove the local existence, c could be chosen 
arbitrary and (A5) gives rjign > (rjit) is continuous). 

For convenience a companion function f " for is defined in the following manner 
(see [m (4.9)]) 

v-{t) = e-^(*-'^V(0) + re~^(*-^)5(7"(s),w!;„(,))cis for t G [a,a]. (34) 

J a 

and v'^ia + 6) = ip{e) for 6 G [-r, 0]. 

It is shown (see [iHl lemmas 4.6 and 4.7)]) that 

Wv^'it) -w''{t)\\< Pmax{6n,Vb{en)} for t e [a - r,a], n = 1,2, P > 0, (35) 
llw" - w"„(^)||c < <5max{£„,r7ft(£:„)} for t e [a,a], n = 1,2, . . . , Q > 0, (36) 
with P,Q > both independent of t and n. 
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Next, [IHl lemma 4.8] shows that z/ there is a function u : [a — r,a] X such that 
u{t) = hm„_j.oo w^{t) uniformly for t E [a — r, cr], then {t, u{t)) G D and n is a solution to 
dlD, on [a, a]. 

Now we proceed the final part of the proof where we use assumption (A5) instead 
of the Lipschitz property of B. First, we need an estimate for \\B(t^,Uti) — B(t'^,Ut2)\\. 
We consider the splitting of the delay mapping F onto continuous and discrete parts 
F = Fc + Fd (according to ( fTOl) ) and use the local Lipschitz property of Fc (due to (A4), 
see (l23l) ). More precisely, using the local Lipschitz property of mapping G ( |32l) . ( I23l) and 
the form ([3]), we get 

\\B{t\uti)-B{t'',Ut2)\\ 
< PR {\t^ - t^\) + Lg,rLf,,r ■ \t^ - + Lg,r{1 + Lf,,r) ■ \ \uti - Ut2\\c 

+ LG,R\\Fd{t\ut^) - Fd{t\ut2)\\. (37) 
Since \t — 7"(t)| < Sn it follows that 

|7"(t) - 7™(t)| as n, m oo uniformly for t G [a,cr]. (38) 

Using ([25D, (ISnD, consider ^ > such that ||w"(t)|| < R for all n > 1 and t G [a — r, cr]. 
Then (|37|) implies 

l|i?(7"(3),<.(.))-i?(7™(3),<™(.))ll 
< z/^ (|7"(5) - l^{s)\) + ■ \r{s) - T{s)\ + Lg,;j(1 + L^„r)| - <™(,)||c 

+LG,i?||F,(7"(s), <„(,)) -Frf(7™(s),w;™(,))|| < [estimates dMD give] 
< La,R{l + L^.,ii)||< - vTWc + Lg,;j||F,(7"(s), <„(,)) - F,(7"^(s), <„(,)) || + e^,™, 
where — )■ as n, m — )■ oo (due to (l38l) ). Using (!34l) . we obtain 

I |t;'^(t) - t^'"(t) 1 1 < f MLa^Ril + L^^ii) I \v: - vTWc ds 

J a 

+ rMLG,;i||F,(7"(s), <„(,)) -F,(7'"(s),<™(,))|Ms + 5.,^, (39) 

^ a 

where in^m — )■ as n, m — >• oo. 

Now our goal is to show that the second integral in fl39l) tends to zero as n, m — oo. 
Let us denote by ip the extension function (^(0 + 6^) = ip{6) for 6' G [— r, 0] and ip{s) = ip{0) 
for s G (a,(j]. We remind that cr — a < T^jg^ (see remark A above). An important 
consequence of (A5) is that Fd{t,Ut) = Fd(t,ipt) for all t G [a, a] and any continuous 
function u : [a — r,(T] X , satisfying Ua = Since all w", by construction, are 
continuous and satisfy = (p, we arrive to the property Fd(t,w^) = Fd(t,ipt) for all 
t G [a, a]. Hence (see the second integral in (139!) ) 

Cris) ^ ||F,(7"(s),«;!;„(,))-F,(7™(s),«;™ (,))|| = ||F,(7"(s), ^,„(,))-F,(7"^(s), ^,™(,))||. 

The continuity of Fa and (!38l) give G^''"(s) — )■ as n, m — )■ 00 for all s G [a,cr]. Since 
(j)J'™(s) is bounded, the classical Lebesgue-Fatou lemma implies that 

t 

G^'"^{s) c?s — )■ as n, m — )• 00. 
13 



The last property gives (see (15^ ) 

\\v-{t)-v^m < f MLa,n{l + LF^,R)\K-vT\\cds + e^,m, (40) 

J a 

where en,m — )■ as n, m — )• oo. 

The rest of the proof follows [18], page 43]. Defining qnmif) = niax{| |t>"(s) — f™(s)|| : 
a — r < s < t} we see that for each t G [a, a] there is an a{t) G [a — r, t] such that 

qn,m{t) = I l^"(«(t)) -v'^{a{t))\\< / MLa,n{^ + Lp^^n) \ |< - v^Wc ds + 

J a 

r-ait) 

< / MLG^ji{l + LF,,R)qn,7n{s)ds + en,m- 
J a 

Gronwall's inequality along with the fact that Sn^m — as n, m — )■ oo shows that qn,m{t) — ^ 
as n, m — )■ oo, and hence {v^{t)}^^i is uniformly Cauchy on [a — r, a]. This implies that 
{w"'{t)}'^^i is uniformly Cauchy on [a — r, a] and hence initial value problem ([T]), ([2]) has a 
mild solution on [a—r, a] (see the discussion above and [iHl lemma 4.8]). The uniqueness of 
solution is due to (A5) and provided by theorem 3. The standard continuation arguments 
give solutions defined on a maximal interval. The proof of theorem 5 is complete. I 

The following important corollary remains valid in the presence of SDDs. 

Corollary (c.f. [LS^ page 18]). Suppose K is a closed, convex subset of X and all the 
assumptions of theorem 5 are satisfied with D{t) = K for all t > a. Suppose further that 

(a) T{t) : K K forty and 

(b) lim;,_o+ ^^(¥^(0) + hB{t, ^y,K)=0 for {t, y;) G V. 

Then ^ has a unique noncontinuable mild solution u on [a, h) for some b > a 
and u{t) G K for all t G [a — r,b). 

Since we are interested in models from biology (see the examples above) the following 
remark is of prime importance for us. 

Remark (c.f. flEl page 7]). Consider system If K = [0, oo)'", then condition 
(a) of the previous corollary holds and condition (b) holds only in case G = {G^)^ is 
quasipositive: if k G {1, . . . , m} and {t, ip) G [a, oo) x C([— r, 0]; C{pt) ) with (p\0, x) > 
for all —r < 9 < 0,x E fl and i = 1, . . . ,m, then y?*(0, ■) = implies G^t, F^t, ip{-, x))) > 
for all X G fl. This condition gives criteria to determine if solutions to ^ remain 
nonnegative if they are nonnegative initially. 
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